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ABSTRACT 


A Nike -Tomahawk sounding rocket was launched near local noon 
on November 8, 1966 from Wallops Island with a triaxial electric 
field experiment using long cylindrical dipoles for sensors. Twenty 

j. 

ELF whistlers with dispersion constants ranging from .15 to .85 sec 2 were 
detected and identified as having propagated up to the rocket through 
the E region from the earth-ionosphere waveguide. The data extend to the 
smallest dispersion ever encountered and the lowest frequencies observed 
on ionospheric whistlers (60 Hz). The waves were eliptically polarized, 
becoming more nearly circular at higher altitudes. Analysis showed 
good agreement between the experimental results and quasi- longitudinal (QL) 
propagation theory for amplitude, dispersion, and polarization characteristics 
of the wave once one is away from the D to E region boundary. No 
appreciable effects of negative ions or protons were observed in the dispersion 


of the whistlers. 



INTRODUCTION 


A long history exists of VLF and ELF measurements from ground based 
instruments, (see Hell iwe 1 1 , 1965). Storey , (1953), postulated that 
whistlers were created by atmospherics and propagated in the circularly 
polarized mode roughly along magnetic field lines in space from the 
northern to southern hemispheres and vice versa. Recently satellite and 
rocket measurements have added to the data and given a better understanding 
of VLF emissions and ion effects on the dispersion of waves propagating in 
the whistler mode. Most previous measurements, including the first by 
( Cain , et. al. , 1961) have been made in the magnetic mode (e.i., Gurnett , 
et. al., 1965; Gurnett , 1966; Smith and Angerami , 1968; Helliwell , 1969 
(review); and the references contained therein). Measurements of the electric 
field are relatively new and have been confined to single or double axis 
measurements ( Gurnett , et al. , 1969; Shawhan and Gurnett , 1968; Scarf , 
et. al. , 1968; Maynard and Heppner , 1969; and the references contained 
therein) , or to fixed frequency propagation experiments ( Storey and 
Cerisier , 1968). 

The rocket experiment reported here represents the first long dipole 
3 axis measurement with wide band telemetry such that the complete wave 
characteristics are determined. Hence, it is possible to check the 
experiment against propagation theory (see Helliwell , 1965) for polarization, 
dispersion, amplitude and direction of the waves. It is also apropos to 
study whistlers in the E region because the angle of injection into the 
E region of the whistler will determine whether it can be guided along 
the field lines, and the small amount of dispersion which occurs prior to 
entry into the ionosphere can be observed. 
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THE EXPERIMENT 

The Nike-Tomahawk sounding rocket (NASA 18.04) was launched from 

Wallops Island on November 8, 1966 at 12:32 hours local time. The 

instrumentation included a three axis electric field experiment using 

the double floating probe technique (A ggson , 1969) in which the difference 

in floating potential between the two probes is monitored. The potential 

difference is the product of the length between the midpoints of the 

— £ > 

probes and the component of the electric field and V x B vectors along 
the axis of measurement. The differential outputs of each axis were 
amplified by 100 with capacitively coupled amplifiers (to separate the 
A. C. components) and used to modulate high frequency, wide-band VCO's 
for wide-band, real time data. The response (see Figure 1) was rolled 
off on the low frequency side by the amplifier with the 3db point at 90Hz 
and on the high frequency side by the response of the ground station 
(10 kHz max. in the data presented in this paper) . 

The cylindrical antennas were 30 feet long with the inner 15 feet 
insulated from the plasma. This gave a 15 meter baseline between the 
mid-points of the active elements. The antenna configuration is shown in 
Figure 2. Thus, the sensitivity was such that 1 mv. input (or 67 pv/meter 
field strength) produced 0.1 volts output from telemetry, and noise levels 
down to about 10 pv/m. could be detected neglecting noise contributions 
from the payload and ground station. The rocket was despun before deployment, 
and the spin was further reduced by deployment. The final motion was 
a slow tumble with a period of about 3 minutes; hence, the spatial 
orientation of each axis varies slowly in time. 

The use of these long dipoles should be contrasted with the short 
base line experiments of Scarf , et.al., (1966) and Scarf , et. al. , (1968). 


Assuming an index of refraction of the order of 100, the wave length of 



-3- 


a 3kHz wave in the ionosphere is 1000 meters which is two orders of 
magnitude greater than the dipole antenna base line. The E = Vmeas./d 
relationship is thus unambiguous. For electrostatic waves ( Scarf , et. al, 
1968) this relation is not necessarily true, because electrostatic waves 
can have wave lengths of the order of or greater than the Debye length 
which is of the order of centimeters in the ionosphere. For electrostatic 
waves of these wavelengths, the electric field is given by a frequency 
dependant relation E = V meag ^ where X is the wavelength (X « d) . The 
only A. C. fields measured, however, were easily recognized as whistler-like 
disturbances from their dispersion characteristics; thus only the whistler 
mode need be considered here. 

THE RESULTS 

The launch time, chosen to optimise the D.C. objectives of the 
experiment, was at midday when minimum whistler activity would be expected. 
During the flight, which lasted approximately five minutes above the 100 km 
level, no whistlers were observed above 100 km with dispersion constants 
long enough to have propagated through the magnetosphere from the 
southern hemisphere. Thus, all of the twenty observed whistlers had to have 
come from the local ionosphere. As the dispersion constants of these 
whistlers ranged from 0.15 to 0.82 sec %, they are all short fractional hop 
whistlers coming from the local earth ionosphere waveguide. 

The noise levels observed from each of the three axes were in general 
less than 50 pv/m, if this noise is interpreted as having wavelengths 
greater than 10 meters. If the source of this noise was electrostatic with 
wavelengths of the order of 10 cm. , this upper limit of background electric 
field strength would correspondingly be less than 5 mv/m. Figure 3 
illustrates three examples of background noise taken from both up and down 
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legs and from near apogee. Part of this noise is from the commutators 
and other electronic equipment on board the vehicle. 

Figure 4 depicts raw data (signals above lKHz attenuated by the 
chart recorder) for two whistlers. The lower whistler observed at 
171.2 km, altitude had a maximum signal amplitude of 450 pv/m resulting 
in a signal-to-noise ratio of 10 to 1. The typical whistler observed 
had a maximum field strength of 100 to 200 pv/m. This whistler had the 
lowest observed frequency components, down to 60 Hz, and was basically 
eliptically polarized. Having just entered the ionosphere (111.7 km), 
the upper whistler had the shortest dispersion of the observed whistlers. 

One also notes that it was almost completely linear in polarization. 

In some cases the rocket was oriented so that one component of the 
whistler was below the noise level. An example of this is shown in 
Figure 5 in which the raw data from the two remaining axes are plotted 
against time. No attempt was made to smooth the data. The high degree 
of circular polarization is evident. The sense is right-handed corresponding 
to the ordinary wave (see Helliwell , 1965). 

In a number of cases, but not in all, frequencies of several kHz 
were observed prior to the ELF portions. The range between 700 to 800 
Hz and 3 kHz was naturally attenuated below the noise level, as explained 
later. The amplitude of the VLF portion sometimes ranged higher than 
800 pv/meter with ratios of the maximum VLF to maximum ELF portions 
being as high as 8 to 1 and typically near 2 to 1. A spectral contour plot 
of a whistler observed at 196 km is shown in Figure 6 to illustrate the 
observed dead band. 

At certain times during the flight modulation from commercial 
broadcast stations was received on the AC channels. The stations would 
fade in and out as the rocket moved through its trajectory. A theory 
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by ( Kawashima , 1968) explains this as rectification by the plasma at 
the plasma frequency. Hence as the vehicle passed through a region 
where the plasma frequency equaled the broadcast frequency, the signal 
would be rectified and the modulation detected by the antennas. 

ANALYSIS 

1. Amplitude 

An averaged amplitude spectrum was obtained using the data from 
four whistlers recorded between 170 and 211 km. having both ELF and VLF 
frequency components. To equalize the contributions each whistler was 
normalized in amplitude to unity at the ELF maximum. The resulting 
spectrum is shown in Figure 7. This spectrum has also been corrected for 
the amplifier response curve of Figure 1. One notes the ELF maximum to 
be at 300 Hz with the amplitude gradually falling off at lower frequencies. 
The aforementioned dead band is seen with the VLF maximum being above 
5 kHz. 

The shape of the spectrum is determined by the source, attenuation 
from propagation in the earth-ionosphere waveguide and attenuation from 
propagation in the ionosphere. It has been shown theoretically and 
by ground based measurements ( Alpert , et. al. , 1966) that the dead band 
between 800 and 3000 Hz is the result of propagation over large distances 
through the earth-ionosphere cavity before entering the ionosphere. This 
is due to attenuation in the waveguide where the cavity spacing is close 
to one half wavelength ( Helliwell , 1965) . Curves for the spectrum of the 
electric field after having propagated various distances from a constant 
intensity source are given by ( Alpert , et. al. , 1966). 

Full wave calculations have been made by (Aksinov and Li shin , 1967) 

for the propagation losses through the lower ionosphere, considering a 
multicomponent plasma between 50 and 200 km, a vertically incident wave 
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and a vertical magnetic field. For comparison here, their resulting 
transmitted energy ratio was converted to a relative electric field 
ratio by taking its square root and multiplying by the fourth root of 
the frequency, to correct for the change in index of refraction. By 
taking the product of this ratio with the curve obtained by ( Alpert , 
et al. , 1966), for a propagation distance of 1000 km, the solid curve 
in Figu;re 6 was derived. Had 2000 km been used, the maximum in the ELF 
region would have broadened and shifted to below 200 Hz and the dead band 
would have been wider. 

Considering the many approximations involved including a constant 
amplitude source, vertical incidence on the ionosphere and vertical 
magnetic field, the agreement with the data is good. From the agreement 
of the data and calculation, it is probable that the whistlers traveled 
over 1000 km before entering the ionosphere. In support of this, there 
was no thunderstorm activity near Wallops Island that day. It is noted 
that the large ratio of VLF and ELF amplitude in the calculated curve 
was seen in some cases although it was not typical. 


2. Dispersion 

Since the frequencies observed here are well below the nose 
frequency the dispersion of the whistler can be described by the relations 

D = r/T (1) 


where T is the group delay, f the frequency, and D is the dispersion 
constant. Hence, the dispersion constant can be shown to be the slope 

j, 

of the T vs. f 2 plot of the measured wave. T may also be related to the 
parameters of the medium in the Storey or quasi longitudinal (QL) approximation 
of Appleton Hartree theory by 


i 

2 c 


fo 

l/f (-f * G )‘ 


ds 


( 2 ) 


T = 


'pprru 
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where f is the electron plasma frequency, is the electron 
gyro- frequency, and 9 is the angle between the wave normal and the 


magnetic field. Thus in the low frequency approximation 

D ~ — ( — ds 


2. C JrnTt? J-f h cos 6 * 


( 3 ) 


This includes only dispersion from electrons and doew not include ion 
effects (see Storey , 1953 and Helliwell , 1965). 

Upon entering the lower edge of the ionosphere the index of 
refraction in the ELF range increases rapidly by two orders of magnitude, 
by Snell's law, a wave intersecting this lower boundary will be refracted 
nearly vertical, hence the wave normal can be assumed to be in the 
vertical direction for all the whistlers above the lower E region. This 
fixes the angle 9 to be the magnetic dip angle which is 20.5° at Wallops 
Island. The ray path or path of energy flow needed for integration of 
Equation 3 will not necessarily follow the wave normal, but as shown in 
the appendix the ray path direction is to a good approximation, at an angle 
from the wave normal whose tangent is 1/2 of tan 9. The integration in 
Equation 3 was carried out numerically starting at 100 km using an inverse 
cube fall-off of magnetic field and the electron density as measured by the 
Wallops Island ionosonde. Since the slope of N g vs h was monotonically 
increasing, the resulting values of N e are unambiguous; however, small 
errors in the altitude assigned are inherent in the measurement. A 
second calculation was done as a check using the electron density measured 
on board the vehicle by the VLF Doppler shift technique of ( Maeda , et al., 
1964) . The results are plotted as the solid and dashed lines respectively 
in Figure 7. 

Since whe whistlers were too short in duration to get an accurate 
measurement of the dispersion from a sonogram type of analysis, they were 
displayed on strip charts and the period of each half cycle of the wave was 
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measured. The time corresponding to this period was assumed to be the 
center time of the half cycle. Zero crossings of the wave were used as 
the limits of the measurement interval. The data for each whistler observed 
were plotted in the form 1/,/f versus t. Three examples of these curves 
taken from different attitudes are shown in Figure 8. The slope of the 
curve was then found using a linear least squares fit; the reciprical of 
this slope is the dispersion constant. 

A systematic error is introduced in the lower E region whistlers 
by the fact that the frequency and amplitude are changing within a half-cycle 
of the wave, hence the assigned time is not the correct time for the 
measured period. Analysis showed that this effect resulted in an error 
of the order of 10% or less on the whistler with a dispersion consistant 
of 0.15 sec \ and correspondingly smaller errors on those with larger 
dispersion constants. The effect is to increase the dispersion constant. 

The resulting value of the dispersion of each whistler is plotted 
in Figure 9. The error bars are a 3cr variation taken by using the 
difference between the slope of the line between each point and the initial 
time and the least squares slope and do not include the systematic error. 

The agreement with the above calculated curve is good down through 130 km 
and deviates slightly below that. This can be accounted for in that: 

(a) there is some dispersion at low frequencies in the earth- ionosphere 
waveguide and in the D region, (b) that the approximations used in the above 
calculations breakdown within a few wavelengths of a sharp boundary, (c) 
and that the systematic error discussed in the previous paragraph exists. 

For eight of the whistlers the VLF portions were analyzed and 
the corresponding dispersion constants were found. The ELF and VLF results 
for these eight whistlers are summarized in Table 1. Variations in the 
dispersions can be expected from the higher frequency of the VLF portion 
(resulting in larger errors in the low frequency approximations) and from the 
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effect of protons. ( Barrington and Nishizaki , 1960) approximated these 
effects by the relation. 

f - = _L f i. | L + >z± 

where fH is the electron gyro frequency 
fiH is the proton gyro frequency 
e is the ratio of the number of 
hydrogen ion to electrons. 

The net result is that a slightly larger dispersion constant would be 

expected at VLF frequencies compared to the ELF frequencies. 

Although the data are scattered, the average ratio of VLF to ELF 
dispersion was found to be 1.01 with a 3a value of .06. If we assume 
a constant magnetic field and that e is an average value over the 
integration region, then the parenthetical factor can be removed from 
the integral. The expected ratio was then calculated for an e of zero 
and found to be 1.005. The effect of the protons would be to slightly increase 
the ratio. ( Brinton , et.al., 1969), have found the number of protons at 
250 km. to be 2 X 1CP /cm 3 and falling off rapidly below that altitude. 
Considering the accuracy of the measurement, and that n is of the order 
of 1 (f /cm 3 e must be assumed to be effectively zero. 

The possibility of a strong interaction of whistlers with negative 
ions in the lower ionosphere was discussed by Shawhan , (1966). Due to the 
resonance of the rotation of circularly polarized whistlers with the 
circular rotation of negative ions at their gyrofrequencies , Shawhan 
showed that the Appleton Hartree dielectric constant diverges at a 
frequency near the negative ion gyrofrequency, if negative ions are present 
in the lower ionosphere. Shawhan chose a rather speculative ion He for 
his calculations; the existance of this ion has not been detected in the 
lower ionosphere. However, negative ion mass spectrograph experiments have 
detected 0“ ions in low quantities in the lower E region (Narcisi , et. al., 



1969) and (Johnson and Heppner , 1956), and it would be of considerable 
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interest to repeat Shawhan's calculations for that ion. It is somewhat 
surprising that these ions do not seem to produce a measurable effect 
below 100 Hz (see Figure 8); however, the response of the amplifier 
prevents looking directly at the 0" gyrofrequency. 


3. Polarization 

The polarization of all the whistlers except those in the 
lower E region was eliptical and in some cases nearly circular. Figure 5 
is illustrative of the whistlers observed near apogee. Since the whistlers 
may have a component in the wave normal direction and magnetic aspect 
only orients the vehicle on a cone, definition of the wave normal direction 
and ray path becomes complicated. 

The polarization equation is given by ( Budden , 1961) as: 


p = - c ^ 

' i-v ? V 


±i. 


% 


4 


2 Y l (,-X-lZ) 


rl 4- j 


where the a direction is that of the wave normal and 

Y t = Y sin 8; Y^ = Y cos 9 


V _ fH 

Y = j 5 * 


X = 


fo 2 


T~ 


z = 


(5) 


9 = angle of B from z axis 

v = collision frequency 
As X > > Y > Z for the observed frequencies and ionospheric conditions 
existing in the E region and above the expression reduces to 


p = + i 


( 6 ) 


or pure circular polarization. There can also exist a component 
of the electric field in the wave normal or z direction given by 

_E a _ -i Yr(^-i) 

( / - X - <- zt ) 

where n is the index of refraction. Considering the parameters at 
210km the index of refraction was calculated to be 373 for a 300 Hs wave. 


Correspondingly, an Es/Ey ratio of about 1.2 would be expected from 
equation 7. 
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The data from two whistlers where one component of the wave 
was zero were plotted to obtain the polarization elipse. The elipse was 
rotated about its minor axis into a circle. The normal to this circle was 
assumed to be the wave normal and was then compared with possible positions 
of the local vertical (as determined by the magnetic aspect) in Figure 10. 
The noted coordinates are those of the vehicle. The agreement is within 
the accuracy of the measurement. The ratio of the component along the 
normal to thatperpendicular was in Figure 10a, 1.3, and in Figure 10b, 

0.7. It is noted that the index of refraction on which the 1.2 estimate 
was made is strongly dependent of frequency and that, as one can see from 
Figure 4, the polarization plot is not precise; hence, the agreement is 
within expected variations. 

CONCLUSIONS : 

The agreement between quasi- longitudinal (QL) propagation theory and 
experiment for the amplitude, polarization and dispersion characteristics 
of waves propagating in the whistler mode has been demonstrated in the 
lower ionosphere. Whereas the simple approximations become crude near 
the D to E region boundary due to the rapidly changing environment, 
the agreement in the upper E region and lower F regions is good. The 
data extends down to the smallest dispersions ever encountered (0.15 sec%) 
and the lowest frequencies observed on atmospheric whistlers (60Hz). 

The agreement of the wave normal data and the dispersion constant data 
with theory provides a check on the ray path behavior in the lower 
ionosphere. No appreciable affects of negative ions was observed in 
the dispersion of the whistlers down to 60 Hz; although changes due to 0" 
might be expected from the work of Shawhan , (1966) . 

The close experimental- theoretical agreement is an indication that the 
experiment can be used to accurately monitor electric fields through 
ELF frequencies . 
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APPENDIX 


Calculation of the ray path angle. 

Let 9 = angle of the wave normal with respect to B. 

a - angle of the ray path with respect to the wave normal. 
Then from Helliwell , 1965, p. 34 
1 do, 


tan o' = 


p No- 


where p is the index of refraction given by Helliwell on p. 29 
for the Q.L. approximation as: 

2 _ i , 

^ f (fH COS 0 - f) 

and where f = the plasma frequency 

fH ■ the electron gyrofrequency 

For ELF frequencies in the E and F regions where f «fn and 

f « f Q , equation A2 may be approximated by: 

f 2 

^ “ f fHCOS 0 

Differentiating this with respect to 9 yields: 
dp ( f o ), (% tan 0) 

ae ~ (f a cos e 

Thus combining A4 with Al: 
tan a = % tan 9 


Al 


A2 


A3 


A4 


A 5 
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FIGURE CAPTIONS 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


1: Frequency response of the broad band A.C. amplifier used on 

each axis, 

2: Configuration of the antennas for the triaxial electric field 

monitor, flight 18.04. 

3: Background noise levels from three different times during the 

flight. 

4: The ELF portions of two whistlers. The upper whistler has the 

shortest observed dispersion constant while the lower whistler 
is typical of those observed at higher altitudes. 

5: Two components of a whistler plotted against time (the third 

component was zero) illustrating the eliptical polarization. 

The two dimensional plot is a compressed view of the 312-227 Hz 
portion of the three dimensional plot. 

6: A spectral contour plot of a whistler observed at 196 km. 

illustrating the characteristic dead band between 1 and 3 kHz. 

7 : Average amplitude spectrum normalized to the ELF maximum 

compared with a theoretically predicted spectrum (see text). 

8: Three plots of Jr vs. elapsed time from which the dispersion 

constant was calculated. One of the shortest whistlers, the 
whistler with the lowest frequency components (also shown in 
Figure 4) and a whistler near apogee are shown. 

9: The dispersion constant variation with altitude compared with 

the calculated variation from Appleton-Hartree theory. The solid 
line represents the calculated value using the ionosonde data 
while the dashed line is that using the Doppler shift data 
(see text). 

10: A polar plot in rocket coordinates showning the locus of the 
normal to circular E (should be nearly vertical from Snell's 
Law) compared with the locus of possible JLocal vertical 
positions derrived from the knowledge of B for two whistlers 
near apogee. 



TABLE 1 


Comparisons of ELF and VLF Whistler Dispersions 


Altitude of 

d„ t „ 

D 

Observations 

ELF 

VLF 

171.2 km. 

.426 sec% 

.581 sec 

178.7 

.503 

.515 

183.2 

.550 

.446 

196.3 

.649 

.623 

197.9 

.686 

.700 

206.6 

.743 

.817 

211.4 

.821 

.656 

211.8 

.783 

.790 



AMPLITUDE GAIN 


soo 



FIGURE I 



FIGURE 2 


179 KM UP LEG 


212 KM NEAR APOGEE 


150 KM DOWN LEG 










*^/VHr"Wv > V'nVV^ 









0 ^s^^“«^^^aV*V'avvy*^ 0 




I 

0 


I 

.01 


I I I I 

.02 .03 .04 .05 

TIME (sec.) 


I 

0 


.01 


I I 
.02 .03 

TIME (sec.) 


.04 


.05 


0 


I 

.01 


I l 
.02 .03 

TIME (sec.) 


I 

.04 


I 

.05 


FIGURE 3 




FIGURE 4 





V z (mv) 



FIGURE 5 


FREQUENCY (kHz) 



illililiip 

B ; j : 


B| 


••-H'i ?: 


V/ ’ < J ‘ 

■ •••• •:•■• •>:'••• -v? <.;->V. • 
V" : ••• ...•: •'• - 

y m 




FIGURE 6 







AVERAGE AMPLITUDE SPECTRUM OF 4 
WHISTLERS OBSERVED BETWEEN 170 AND 
21! KM. 18.04 

SOLID CURVE COMBINES CALCULATIONS 
OF ALPERTj El AL, AND AKSENOV AND 




ELAPSED TIME (MSEC) 


FIGURE 8 



HEIGHT (KM) 



DISPERSION CONSTANT (Sec* 


FIGURE 9 




